Users manual for coordinate generation code CRDSRA

Generation of a viable coordinate system represents an important component of an isolated airfoil Navier-Stokes calculation. The manual describes a computer code for generation of such a coordinate system. The coordinate system is a general nonorthogonal one in which high resolution normal to the airfoil is obtained in the vicinity of the airfoil surface, and high resolution along the airfoil surface is obtained in the vicinity of the airfoil leading edge. The method of generation is a constructive technique which leads to a C type coordinate grid. The method of construction as well as input and output definitions are contained herein. The computer code itself as well as a sample output is being submitted to COSMIC

On the application of a hairpin vortex model of wall turbulence to trailing edge noise prediction

The goal is to develop a technique via a hairpin vortex model of the turbulent boundary layer, which would lead to the estimation of the aerodynamic input for use in trailing edge noise prediction theories. The work described represents an initial step in reaching this goal. The hairpin vortex is considered as the underlying structure of the wall turbulence and the turbulent boundary layer is viewed as an ensemble of typical hairpin vortices of different sizes. A synthesis technique is examined which links the mean flow and various turbulence quantities via these typical vortices. The distribution of turbulence quantities among vortices of different scales follows directly from the probability distribution needed to give the measured mean flow vorticity. The main features of individual representative hairpin vortices are discussed in detail and a preliminary assessment of the synthesis approach is made

A compressible solution of the Navier-Stokes equations for turbulent flow about an airfoil

A compressible time dependent solution of the Navier-Stokes equations including a transition turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme. A nonorthogonal body-fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional, and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles of Reynolds number up to one million and low subsonic Mach numbers

Three-dimensional unsteady viscous flow analysis over airfoil sections

A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations

Numerical solutions of Navier-Stokes equations for compressible turbulent two/three dimensional flows in terminal shock region of an inlet/diffuser

The multidimensional ensemble averaged compressible time dependent Navier Stokes equations in conjunction with mixing length turbulence model and shock capturing technique were used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations were carried out: (1) 2 D, steady, subsonic; (2) 2 D, steady, transonic with normal shock; (3) 2 D, steady, supersonic with terminal shock; (4) 2 D, transient process of shock development and (5) 3 D, steady, transonic with normal shock. The numerical results obtained for the 2 D and 3 D transonic shocked flows were compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data

Parametric study of relaminarization of turbulent boundary layers on nozzle walls

By means of comparisons between theoretical predictions and experimental data, the accuracy of a boundary procedure to predict the effect of large streamwise accelerations upon initially turbulent boundary layers is assessed. The boundary layer procedure is based upon simultaneous solution of the boundary layer partial differential equations and the integral turbulence kinetic energy equation. The results of the present investigation show the ability of the procedure to accurately predict properties of boundary layers subjected to large streamwise accelerations. The procedure is used to conduct a parametric study of the effect of free stream turbulence, heat transfer, Reynolds number, acceleration, and Mach number on boundary layers in supersonic nozzles to assist in the design of a quiet tunnel. Results of the investigation show that, even in the presence of moderate free-stream turbulence levels, the boundary layer in the approach section of the quiet tunnel nozzle relaminarizes and becomes thin enough to be removed by a small slot in the nozzle wall. Furthermore, the calculations indicate that it should be possible to maintain a laminar boundary layer for the entire length of the supersonic portion of the quiet tunnel nozzle

Analysis of strong-interaction dynamic stall for laminar flow on airfoils

A compressible Navier-Stokes solution procedure is applied to the flow about an isolated airfoil. Two major problem areas were investigated. The first area is that of developing a coordinate system and an initial step in this direction has been taken. An airfoil coordinate system obtained from specification of discrete data points developed and the heat conduction equation has been solved in this system. Efforts required to allow the Navier-Stokes equations to be solved in this system are discussed. The second problem area is that of obtaining flow field solutions. Solutions for the flow about a circular cylinder and an isolated airfoil are presented. In the former case, the prediction is shown to be in good agreement with data

User's manual for airfoil flow field computer code SRAIR

A two dimensional unsteady Navier-Stokes calculation procedure with specific application to the isolated airfoil problem is presented. The procedure solves the full, ensemble averaged Navier-Stokes equations with turbulence represented by a mixing length model. The equations are solved in a general nonorthogonal coordinate system which is obtained via an external source. Specific Cartesian locations of grid points are required as input for this code. The method of solution is based upon the Briley-McDonald LBI procedure. The manual discusses the analysis, flow of the program, control steam, input and output

Calculation of steady and unsteady airfoil flow fields via the Navier-Stokes equations

A compressible time-dependent procedure for the two-dimensional ensemble averaged Navier-Stokes equations has been applied to the isolated airfoil problem in steady and unsteady flows. The procedure solves the governing equations via the linearized block implicit technique. Turbulence is modeled either via a mixing length or turbulence energy approach. The equations are solved in general non-orthogonal form with no-slip boundary conditions applied at the airfoil surface. Results are presented for airfoils at constant incidence, an airfoil in ramp motion and an airfoil oscillating through a dynamic stall loop. In general, steady converged solutions are obtained within 70 time steps over the range of Mach numbers considered. Comparisons with measured data show good agreement between computation and measurement

A three-dimensional viscous flow analysis for the helicopter tip vortex generation problem

The tip vortex flow field occurring in the vicinity of the tip region of a a helicopter rotor blade is a very complicated three-dimensional, viscous flow phenomenon. The details of the flow in the tip region can have a major effect in determining the generated rotor noise and can significantly affect the performance and dynamic loading of the rotor blade. The three-dimensional viscous subsonic tip vortex generation processes is investigated by a numerical procedure which allows spatial forward-marching integration, utilizing flow approximations from the velocity-decomposition approach of Briley and McDonald. The approach has been applied to compute the laminar and turbulent tip vortex flows for a constant thickness slab airfoil with a square tip, a constant thickness slab airfoil with a half round tip and a NACA 0012 airfoil with a half round tip. The basic mechanism of the tip vortex generation process as well as the prediction of vortex appearance, strength and secondary flow shown by the calculations are in qualitative agreement with experimental results